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ABSTRACT 

We combine recently computed models of stellar evolution using a new treatment of rotation with a 
Bayesian statistical framework to constrain the ages and other properties of early-type stars. We find 
good agreement for early-type stars and clusters with known young ages, including /3 Pictoris, the 
Pleiades, and the Ursa Majoris Moving Group. However, we derive a substantially older age for the 
Hyades open cluster (750 ± 100 Myr compared to 625 ± 50 Myr). This older age results from both the 
increase in main-sequence lifetime with stellar rotation and from the fact that rotating models near the 
main-sequence turnoff are more luminous, overlapping with slightly more massive (and shorter-lived) 
nonrotating ones. Our method uses a large grid of nonrotating models to interpolate between a much 
sparser rotating grid, and also includes a detailed calculation of synthetic magnitudes as a function of 
orientation. We provide a web interface at www.bayesianstellarparameters. info, where the results 
of our analysis may be downloaded for individual early-type (B — V < 0.25) Hipparcos stars. The 
web interface accepts user-supplied parameters for a Gaussian metallicity prior and returns posterior 
probability distributions on mass, age, and orientation. 

Subject headings: 


1. INTRODUCTION 

Stellar ages are important in many areas of astro¬ 
physics, from understanding the evolution of galaxies, 
to chemical enrichment, to the properties of planetary 
systems. They are also usually inaccessible by direct 
measuremen t and must be obtained by fitting to models 


(Soderblom 


2010). The one exception is the Sun, radio- 


metrically "dated from inclusions in meteorites (Connelly 
et al. 2012|. Stars older than ^100 M yr may be dated 
by fitting models of ste llar evolution (Perryman et al. 
1998 Takeda et al. 20071 or by using calibrated empirical 


methods like lithium depletion (Sestito fe Randich 


20051 


or the decline of activity and rotation (Noyes et al.|1984] 


Barnes||2007 Mamajek & Hillenbrand 


20U8| ). These em- 


pirical methods are calibrated at ages of several hundred 
Myr using clusters dated by fitting model isochrones. 

Low-mass Held stars may be reliably dated by the de¬ 
cline of their rotation and activity. FGK stars develop 
an outer convective zone, generating a magnetic field 
and driving a wind that carries awa y angular momen¬ 


tum (Parker 1955 


and become less c 


with age (Wilson 

viding a clock tha t can give ages to 


Glatzmaier 1985). They spin down 
rromosphcrical ly and coronally active 


1963 Hempelmann et al. 1995), pro- 
20% (Mamajek &; 


Hillenbrand 20081. Isochrone fits to coeval clusters are 
used to calibrate the empirical rotation and activity rela¬ 
tions at ages of several hundred Myr. The accuracy and 
precision of gyrochronology ages thus rely on the those 
of isochrone fits. 

Stars more massive than ~1.3 Mq have radiative en¬ 
velopes and convective cores, and retain thei r angu- 
lar momentum throughout the main sequence (Barnes 
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2003). Even individual early-type sta rs are t ypically 
dated by fitting model is ochrones (e.g. Zorec fe Royer| 


2012 Nielsen et al. 2013). Given an initial mass and 
metallicity, stellar evolution models predict the stellar 
radius and luminosity as a function of time. Together 
with a model atmosphere, theoretical isochrones can be 
used to generate a grid of luminosities in each band as 
a function of initial mass, metallicity, and age. A star’s 
distance and fluxes can then be inverted to give a mass, 
metallicity, and age. Unfortunately, the inversion is often 
not unique, with multiple paths in the color-m agnitude 
diagram crossing one another (SoderblompOlO). 

Models of stellar evolution are subject to a very wide 
range of uncertainties, from the treatment of convec¬ 
tion and overshooting, to the composition and opaci¬ 
ties, to the treatment of rotation. Rotation, in par¬ 
ticular, can introduce moderate changes in luminosity, 
color and brightness variations with viewing angle, and 
large, ~25% increases in the main-sequence lifetime as 
rotationall y-induced mixing supplies the core with fresh 
hydrogen (Meynet & Maeder|2000| Ekstrom et al.|2012|). 
The dependence of main-sequence lifetime on rotation 
can broaden the main-sequence and especially the main - 


sequence turnoff (Georgy et al. 2014; Li et al. 


2014), 


mimicking the effect of a spread of ages. The effects of 
rotation are most important for stars that do not brake 
magnetically, i.e., stars >1.3 Mq. These are precisely the 
stars for which isochrone dating is the most important, 
and often the only, way of estimating their age. 

In this paper, we perform isochrone dating of early- 
type stars using a recent set o f stellar mode ls with a shel- 
lular treatment of rotation (Georgy et al. 2013). These 
models cover a range of metaliicities from Z Q = 0.002 
to Z = 0.014, initial rotation rates from 0 to 95% of 
breakup, and masses from 1.7 Mq to 15 Mq. We use 
a much finer grid of nonrotating stellar models (Girardi 


et al. 2002) to interpolate between the rotating models, 
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and to extrapolate to 1.45M 0 and to Z = 0.04. The 
extrapolations rely on the fact that the additional effects 
of rotation in these models depend only weakly on mass 


and metalli city. We use a Roche model of the star (Es- 
pinosa Lara fc Rieutord|2011[) together with the ATLAS9 


ft: 


p(M, Z , r, ft, i) oc p(M)p(Z)p(m)p(i)p(Q)£(M, Z, r, ft, i) . 

( 3 ) 

The appropriate priors on Mo, r, m, and i are all well- 
established: uniform in volume and time, and propor¬ 
tional to the initial mass function (IMF). The prior on 
parallax is thus dp/dm oc m~ 4 , while dp/di = sin i, where 
i = 0 denotes a polar viewing angle. We adopt a Salpeter 
IMF, dp/dM oc M~ 2 35 , appropriate for st ars >1 M 0 
( Salpeter|l955 Kroupa|2001 Chabrier||2003 ). 

An appropriate prior on Z in the absence of other in¬ 
formation about the star may be taken from the metallic- 
ities of nearby young stars. There is some disagreement 
on the chemical composition of the young Solar neigh¬ 
borhood. Th e large spectroscopic sample of Casagrande 
et al. (2011) implies a slightly sub-Solar mean metallic¬ 
ity tor young stars, and a dispersion ~0.1 dex. This may 


disag ree with Przybilla et al. (2008) and Nieva & Przy- 
billa (2012), who compared high resolution, high signal- 


to-noise ratio spectra of twen ty nearby early B s t ars to 
detailed atmospheric models. Nieva & Przybilla (20121 
found abundances almost identical to those of the Solar 


photosphere ( Asplund et ah|2 009), with a scatter of just 
~10%, or 0.03 dex, even when including measurement 
uncertainties. These results are in much better agree¬ 
ment with models of Galactic chemical enrichment and 


model atmospheres (Gastelli & Kuru cz~p 004) to compute 
synthetic photometry for each stellar model as a function 
of orientation. We then adopt a Bayesian framework to 
constrain stellar parameters. 

We organize our paper as follows. In Section [2]we 
describe our Bayesian statistical formalism, while inSec- 
tion[3]we describe the stellar models that we use and the 
effects of adding rotation. In Section [4] we describe our 
use of the Roche model of a star together with model 
atmospheres to compute synthetic photometry. In Sec¬ 
tion [5] we apply our method to four nearby early-type 
stars and to three coeval associations, the Pleiades, Ursa 
Majoris, and the Hyades. We conclude with Section [6] 

2. STATISTICAL FRAMEWORK 

The likelihood of a stellar evolution model is the proba¬ 
bility of measuring the observed fluxes and distance given 
the predicted luminosities in each band: 

£(M, Z,r,m) = p({m j ,m ohs }\{M j (M,Z,T),m}) , (1) 
where M is the model’s mass, Z its metallicity and r 

its age, m is the parallax, rrij are the apparent mag- „ _ 

nitudes in bands j. and Mj are the (absolute) model ( |2U11 ) also show some trends that are unlikely to be real 
magnitudes. The addition of rotation in stellar models 
enables us to include the measured projected rotational 
velocity, usin*, in the likelihood function. Stellar rota¬ 
tion also modifies the main sequence lifetime, apparent 
luminosity, and colors, which all become functions of the 
angular velocity ft and inclination i (i = 0 denoting pole- 
on). These become additional parameters in Equation 
0 upon which the Mj depend. 

Fbr convenience, and for lack of an obvious alternative, 
we assume that the errors in observed magnitudes are 
Gaussian. We neglect errors in the synthetic photometry 
apart from a possible systematic error that we add in 
quadrature to the photometric errors. Such a systematic 
error is best interpreted as a color uncertainty. The error 
in the Hipparcos parallax (which we also assume to be 
Gaus sian ) is ty pically a few percent for nearby stars (van 
Leeuwen 2007), which is equivalent to much more than 
a 1% uncertainty in flux. With these assumptions and 
simplifications, the likelihood function becomes 


mixing in the interste llar medium (R oy fc Kunth 1995 
Chiappini et al. 2003). The data m Gasagrande et al. 


including a correlation of metallicity with stellar mass at 
fixed age (their Figure 16). We suggest a Gaussian prior 
centered on [Fe/H] = 0 with a dispersion of 0.1 dex as a 
very conservative (perhaps overly conservative) choice. 

Our suggested metallicity prior is problematic for 
chemically peculiar early-type stars. These stars are ex¬ 
pected to have unremarkable bulk compositions, but at¬ 
mospheres heavily influenced by diffusi on and/or mag¬ 


netic fields (Preston 1974; Smith 1996). Individual el- 


—2 ln£ = 

bands j 


(Mnod, ^ logiQ ^obs, j) 


(m - m ohs ) 2 (R e q ft sin i - v ohs ) 


( 2 ) 


where M mo( j 0 is the model absolute magnitude in band 
j , m is the parallax in milliarcseconds, R e q ft is the model 
equatorial velocity, and r> 0 b s is the observed projected 
rotational velocity. 

In a Bayesian framework, the posterior probability of 
the model is the product of the likelihood and the prior 
probability of the model parameters: Z, Mo, t, m, i and 


ements and groups of elements can be over- or under- 
abundant by factors of tens or hundreds. Chemically 
peculiar stars are relatively common on the upper main 
sequence, and our results for them must be interpreted 
with caution. Such stars should be perhaps be fitted 
by a stellar model with different bulk and atmospheric 
compositions. 

The prior on stellar rotation is more difficult to deter¬ 
mine. Stellar evolution codes calculate a star’s structure 
throughout its main-sequence life, which can entail a sig¬ 
nificant redistribution of its angular momentum. The 
only free parameter in this case is a star’s initial supply 
of angular momentum, which need not be represented 
by its observed equatorial velocity (and which is shed in 
a magnetized wind by stars <1.3 M 0 ). An appropriate 
prior probability on rotation would be the distribution 
of angular momentum (or equatorial velocity, often the 
only observable) for young early-type stars, which do not 
s hed angular momentum . 

Zorec fc Royer|(2012) have estimated the evolution of 


equatorial velocity tor stars in various mass bins at dif¬ 
ferent stages of evolution. For ~ 1.5-3 M 0 stars early in 


their main-sequence lifetimes, [Zorec fc Royer find peaks 
in the rotational velocities consistent with a bimodal dis¬ 
tribution. The first peak is at low velocities while the 
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second, which accounts for most of the stars, may be 
approximated by a lagged Maxwellian with a mode of 
ft/ft cr it ~ 0.5. Their analysis also hints that fl/fi cr it 
may increase somewhat w ith initial st ellar mass. 

We use the models of Georgy et al. f2013) in our anal¬ 
ysis, which have equatorial Sl/S2 cr ; t somewhat lower than 
the initial values at which the models are tabulated. We 
therefore use distributions in fi/f2 cr it pe aking at some¬ 
what l arger values than those derived by Zorec &; Royer| 
(2012). For simplicity, we use a simple Maxwellian dis- 
tribution with a mode of 0.5 truncated at 0.95. This cap¬ 
tures the essential elements of the data jZorec fe Royer) 
analyze, apart from a few (~0-20%, depending on mass) 
of stars that appear to rotate slowly, and reproduces the 
approximate upper limits to observed values of usini. 

3. STELLAR EVOLUTION MODELS WITH 
ROTATION 

A wide range of stellar evolution models is now avail- 


et al. 2001 


Girardi et al. 2002:j 

Pietrinferni et al. 

20(Hf 

Dotter et a 

1. 2008p. Some models 

, going back more 

han a 


The addition of rotation turns a ID stellar structure into 
at least a 2D structure. This is simplified in practice with 


assumptions like that of “sliellular” rotation (Zahn 1992), 
in which the angular velocity is constant along isobars. 
Rotation is relatively unimportant for stars below ~1.3 
Mq , which spin down rapidly. However, it can introduce 


l arge effects fo r stars with radiative envelopes (Meynet & 
Maeder|20 00), perhaps most notably an increase in time 
spent on the main sequence. 

Stellar modeling includes a very wide range of tun¬ 
able physical parameters governing everything from the 
chemical abundances and their opacities, to the treat¬ 
ment of convection, to the treatment (or lack) of rotation, 
to the model atmospheres. These parameters and treat¬ 
ments vary from one set of stellar models to another and 
can lead to si gnificant variations in colors, luminosities, 
and lifetimes (Lebreton et al.||2014 ). For the rest of t his 
analysis, we adopt the models of Georgy et al. (2013) as 
our primary source, compute synthetic photometry our- 
selves, and neglect uncertainties in parameters like the 
helium fraction or convective overshooting._ 

The rec ent stellar evolution models of lEkstrom et aid 
(2012) and Georgy et al. (2013) include both rotating and 
non-rotating stars, in the former case covering a range of 
initial rotation rates and metallicities. These grids are, 
however, very coarse, and only extend to Z = Zq. We 
therefore use much finer grids of nonrotating models to 
interpolate between the rotating models, and to extrapo¬ 
late them t o higher metallicitie s. We adopt the PARSEC 
isochrones (Girardi et al. 2002 1 for this interpolat ion. Be¬ 
cause we anchor them to the|Georgy et al.|(|2013|) models, 
the details of the PARSEC calculations are unimportant 
here. Any effects that are linearly dependent on mass 
and metallicity will disappear entirely. 

The validity of our interpolation and extrapolation re¬ 
lies on the fact that the modifications to stellar evolution 
provided by rotation depend only weakly on other stel¬ 
lar parameters. In this analysis, we restrict the models 
to M > 1.5 Mq, for which magnetic braking is inefficient 
and rotation is important throughout the main sequence. 
We apply our upper mass cutoff at 11 Mq, which has a 


main sequence lifetime of just a few tens of Myr. 

The rapidly rotating stellar models have main- 
sequence lifetimes that depend strongly on f2, the ini¬ 
tial angular momentum. Interpolating with nonrotating 
models requires us to introduce a dimensionless time, 
2/2ms> the ratio of the stellar age to its main sequence 
lifetime. The parameter 2 ms is a function of stellar rota¬ 
tion: 

^MS = tMS,nr/3 (Af 0) D) . (4) 

We use the PARSEC models to interpolate 2ms, nr be¬ 
tween masses in the Geneva models, and separately inter¬ 
polate /3, which depends very weakly on mass and metal¬ 
licity. We note that the nonrotating PARSEC isochrones 
do differ significantly from the nonrotating Geneva mod¬ 
els. For example, 1.5-5 Mq PARSEC models spend 
~10% longer on the main sequence than the correspond¬ 
ing Geneva models we use here. 

The other effects of rotation include an increase in lu¬ 
minosity and a flattening of the stellar surface, which 
make the apparent luminosity of the star a function of 
viewing angle. We describe the computation of synthetic 
photometry in the following section. To interpolate be¬ 
tween models, we separate the photometry into nonrotat¬ 
ing magnitudes and a rotating correction as a function of 
f2, photometric band, and inclination. These two compo¬ 
nents are additive in units of magnitude (multiplicative 
in units of flux). We then use the PARSEC models to in¬ 
terpolate the nonrotating synthetic photometry between 
Geneva models. We separately interpolate or extrapolate 
the rotation correction term, which, like (3 from Equation 
0 , depends only weakly on stellar mass and metallicity. 
This interpolation in linear in units of magnitude. 


4. COMPUTING THE EFFECT OF ORIENTATION 

Stellar oblateness introduces two effects: a viewer 
along the pole sees a star both larger in area and hot¬ 
ter. In all bands, the star will therefore appear more 
luminous with decreasing inclination i (i = 0 denoting 
a line-of-sight parallel to the polar axis), except when 
i approaches 0° and the gravity darkening effect is se¬ 
vere. The variation of stellar surface brightness with lat¬ 
itude isj^consequence of the von Zeipel theorem (von 


Zeipel 1924), which states that the local effective temper- 
ature depends on the local effective gravity as a power- 
law T e g oc with /3 known as the gravity darkening 
coefficient. The original von Zeipel law states (3 = 1/4 
for stars with fully radiative envelopes. Generally, /3 is 
smaller than 1/4, as discovered by recent interfer ometr ic 


berg et al. 

2006| van Belle et al. 

2006 

)■ 

Althoug 

i comprehensive calcu 

ations of f3 were carried 


out by Claret (1998) for stars at different evolutionary 
stages, we choose t o adop t a two -dimensional model de¬ 


veloped by Espinosa Lara & Rieutord (2011) (hereafter 
LR11), in which the effective temperature profile only de¬ 
pends on the ratio of the equatorial velocity and the Ke- 
plerian velocity w — n/u C rit- This parameter is supplied 


by the stellar evolution models (Ekstrdm et al. 2012). 
With the total luminosity, oblateness, and radius from 
the stellar models, we can compute T e ff and g e g over the 
entire stellar surface. Finally, we use the specific intensi- 
ties from the ATLAS9 models ( Castelli fc Kurucz||2004 | 


to integrate the apparent flux through a series of filters 
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as a function of viewing angle. 


4.1. Computing Fluxes 

We follow Section 2 of LR11 to compute the effective 
temperature map over the stellar surface. The assump¬ 
tion of the model can be simply stated as the following: 
the mass distribution of the star may be described as 
a Roche model, and the energy flux is a divergence-free 
vector that is almost anti-parallel to the effective gravity. 
The effective temperature profile at a particular polar 
angle 9 can then be expressed as 


T eS = [ ~^g + w 4 r 2 sin 2 9 - 

1/4 


2w 2 sin 2 9\ 1 ^ 8 


V 4 7rcrR2 q 


F 1 / 4 


(5) 


in which r is the radius at the stellar surface in units of 
the equatorial radius R eq , w is the ratio between the 
equatorial velocity and the Keplerian velocity v/v CI g, 
and F w is a correction factor that takes into account 
the difference between the vector angle of energy flux 
and the effective gravity. When the rotation is slow 
(-R P oie/-Req > 0.95), F v , = 1, and this model reduces 
to the von Zeipel case. 

We solve for F w numerically using a combined Newton- 
Raphson and bisection method following Equations (24)- 
(28) of LR11. The stellar surface is assumed to be an 
oblate ellipsoid with a = R eq and c = R po i- The pro¬ 
jection of the stellar surface onto the plane of the sky 
with an inclination i is an ellipse with a = R eq and 
b = (R 2 q R 2 ol )/(i? 2 q cos 2 i + i? 2 ol sin 2 i). 


We use the ATLAS9 stellar atmosphere models 


(Castelli & Kurucz 
sities, and the Tyc. 


2004 ) to co mpute the specific intcn- 
10 , 2MASS (iCutri et al.l 20031), and 


Johnson-Cousins filter response curves to compute pho¬ 
ton fluxes at a distance of 10 pc. We adopt the r ecently 
revised Tycho responses (Bessell & Murphy[2012|), which 
differ slightly (b y <5 mmag in Ijt — Vt ) from those origi¬ 
nally published '(ESA 1997). T he a bund ances in t he AT - 
LAS9 models are relative to Anders & Grevesse (1989), 
which gives a Solar metallicity Z = 0.019. This is higher 


by_0.10-0.13 dex than the JSL. values adopt ed by |Gi- 
rardi et al. (2002) and Ekstrom et al. ([2012 ) based on 


the abundances measured by Asplund et al.| (2009). We 
therefore use the ATLAS9 models with [Fe/H] = — 0.1 
for the isochrone models at solar metallicity. 

Our adoption of Zq ss 0.014 could introduce system¬ 
atic effects when using spectroscopic metallicities that 
assume a different ^0- As a result, [Fe/H] = 0.1 in this 
work could correspond in some wa ys to a spectroscopic 


[Fe/H ] = 0.0. Our scaling of the Anders & Grevesse 
(1989) photospheric composition does, however, mean 
that the composition of the stellar atmosphere roughly 
matches the bulk composition of the star. 


4.2. A Parametrized Fit 


of Georgy et al. (2013) and PARSEC models at the same 
mass, metall icity, and fract i on of main sequence lifetime. 
Some of the Georgy et al. (2013) models are computed 
with zero initial angular momentum. These models are 
not identical to the corresponding PARSEC models; for 
example, nonrotating PARSEC ~2 M 0 stars typically 
spend ~10% longer on the main sequence than the cor¬ 
responding nonrotating Geneva models. Even at zero 
rotation, the difference between Geneva and PARSEC 
model magnitudes can be significant. 

We compute the magnitude corrections as a function 
of viewing angle p = cos i (p = 1 corresponding to pole- 
on). We then fit a polynomial in /i for each band and 
rotating model, 


M(p = cosi) = Mparsec + ^2, ai ^ • (6) 

»=o 

Figure [l] shows the residuals from this fitting function for 
the Johnson-Cousins U passband as a function of T e g and 
obliquity at an equatorial log g = 4. We obtain excellent 
agreement with n = 5, with residuals <1 mmag for all 
but the coolest models at the fastest rotations. These 
residuals decrease towards longer-wavelength passbands. 
With n = 5, they are typically at least an order of magni¬ 
tude smaller than the photometric measurement errors, 
and likely smaller than errors resulting from interpola¬ 
tion between isochrones. 

The ai in Equation © are functions of the band, stel¬ 
lar mass, rotation rate, metallicity, and age, and are log¬ 
arithmic in the stellar flux. We linearly interpolate (and 
even extrapolate, in metallicity), these parameters onto 
the fine nonrotating PARSEC grids. In doing so, we take 
advantage of the fact that while the stellar models them¬ 
selves may depend strongly on physical parameters like 
mass and metallicity, the at correction coefficients have 
a much weaker dependence. 

With the effect of rotation reduced to the parameters 
of Equation ([6| , we can compute the effects of both rota¬ 
tion and orientation on observed colors and magnitudes. 
Figure [2] shows both nonrotating and rotating isochrones 
(with tr/Cl clit = 0.5) at a series of representative ages. 
A rotating model at fixed age actually covers a line in 
the color-magnitude diagram, as much as ~0.1-0.2 mag- 
ntiudes in V and 0.05 magnitudes in B — V (the navy 
blue lines in the figure). At very young ages, the effect 
of rotation is modest. At older ages, however, a rotat¬ 
ing stellar model will overlap a significantly younger and 
slightly more massive model in color-magnitude space. 
The combination of an increased luminosity at a simi¬ 
lar B — V color with a longer time spent on the main 
sequence can have a very large effect on isochrone ages 
using main-sequence turnoff stars. This is especially ap¬ 
parent for the Hyades open cluster, which we discuss at 
the end of the next section. 

5. RESULTS 


Our method allows us to compute the magnitudes in 
all bands as a function of orientation. We wish to use 


the fi ne grid of nonrotating PARSEC isochrones (Girardi 


et al. 2002) to interpolate within the coarse grid of ro¬ 


tating Geneva models. We therefore compute the differ¬ 
ence in magnitude between the rapidly rotating models 


In the following subsections we present the results of 
our isochrone analysis for a selection of A and B stars. 
We apply it first to a series of three individual stars of 
known age and to one star, k And, with an imaged sub- 
stellar companion and a controversial age. We have gen¬ 
erally excluded close visual and spectroscopic binaries 
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Oblateness Oblateness Oblateness Oblateness 



Fig. 1.— Distribution of the residuals from Equation § for the Johnson-Cousins /7-band. With a fitting polynomial of order n = 5, the 
root-mean-square residuals are <1 mmag for all but the lowest temperature, most rapidly rotating stars. This is much less than pole-equator 
differences reaching 0.1 mag or more, and less than errors due to interpolation and measurement errors. The residuals are smaller still 
(i.e. the fit is better) at the longer-wavelength passbands that we use. 



Fig. 2. — Isochrones for three stellar masses, with (blue) and 
without (red) rotation at fixed f2/f2 cr it = 0.5; several ages are 
noted in Myr. The rotating models have been plotted only at rep¬ 
resentative ages to show the effects of orientation. The lines show 
the positions in color-magnitude space as a function of orientation. 
At older ages, a rotating model will overlap a significantly younger, 
slightly more massive nonrotating star in color-magnitude space. 


but have not excluded any stars based on the perfor¬ 
mance of our fitting. We then apply our isochrone-based 
analysis to three stellar clusters with previously inferred 
isochrone ages, the Pleiades, the Ursa Majoris moving 
group, and the Hyades open cluster. 

We restrict our analyses to the Tycho-2 Bt and Vr 


photometry (Hpg et al. 2000). These data have been 
well-calibrated and offer precisions of ~0.01 mag down 
to Vr ~ 8 over the entire sky. 2MASS JHK S photometry 
could also be used, but would require a careful analysis 
demonstrating its calibration relative to Tycho at levels 
of <0.01 mag, and a demonstration that the model at¬ 
mospheres are adequate over such a wide range of wave¬ 
lengths. We defer such an investigation to a f uture paper. 


We u se Hipparcos parallaxes throughout (van Leeuwen 


20071, an d use usin i measurem ents collected and cali¬ 
brated by |Zorec & Royer| ([20121. 

In all of our fitting, we apply floors on the measure¬ 
ment errors of all parameters, which we add in quadra¬ 
ture with the actual reported uncertainties. These error 
floors are 0.005 magnitudes in all bands (^0.5% in flux), 
which are best thought of as color errors, and 30 kms -1 
in usini. They are designed to minimize the impact of 
the coarseness of our model grid and to account for some 
uncertainties in the synthetic photometry. The Geneva 
models are only computed at nine initial angular mo¬ 
menta; the spacing between them corresponds to ~20- 
50 kms -1 in equatorial velocity. The spacing between 
neighboring models in luminosity is typically ~2% (and 
typically much less in color). The difference in synthetic 
colors between the original and recalibrated Tycho pho¬ 
tometric system produces a systematic shift in Bt — Vr 
of ^0.004 mag, while errors in parallax are >2%. 

When fitting an ensemble of stars that share an age and 
chemical composition, we multiply the two-dimensional 
posterior probability distributions for each individual 
star. This relies on the stellar models providing ade¬ 
quate fits to all stars, and on these fits being consistent. 
A rigorous statistical statement of consistency depends 
on the actual distributions. Roughly speaking, however, 
consistency requires that the ler contour of the multi¬ 
plied posterior distribution be at a level comparable to 
the product of the ler contours in each individual distri- 
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Fig. 3.— Posterior age probability distributions for four Hippar- 
cos stars, each with three metallicity priors: [Fe/H] = 0.1 zb 0.05 
(red histograms), [Fe/H] = 0 =b 0.05 (blue histograms), and 
[Fe/H] = —0.1 d= 0.05 (green histograms). HIP 13209 (third panel) 
has a binary companion of unknown spectral type which we neglect 
in the isochrone fit. The shaded light blue regions show the moving 
group ages of the top three stars. Slightly super-Solar metallicities 
allow an excellent match in all cases. HIP 116805’s age is contro¬ 
versial; our analysis suggests that youth cannot be ruled out. 


bution. All of our cluster probability distributions pre¬ 
sented below easily pass this test. 

5.1. Single Stars 

We first apply our isochrone-based analysis to four in¬ 
dividual stars:'HIP 27321 (/3 Pic), 13209, 115738, and 
116805 (ft And). The first three of these stars are con¬ 
sensus members of the coeval moving groups /3 Pictoris 
(HIP 27321) and AB Doradus (HIP 13209 and 115738). 
HIP 116805 is a proposed member of the Columba mov¬ 
ing group, but this identification, and the star’s age, are 


controversial. Figure [3] shows the age posterior probabil¬ 
ity distributions for each of these stars for three Gaus¬ 
sian metallicity distributions: [Fe/H] = 0.1 ± 0.05 (red 
histograms), [Fe/H] = 0 ± 0.05 (blue histograms), and 
[Fe/H] = —0.1 ± 0.05 (green histograms). The shaded 
blue regions show the moving group ages. We discuss 
each star in turn below. 

HIP 27321: The nearby A star /3 Pictoris h osts a de¬ 


bris disk and a low-mass substellar companion (Lagrange 


toris moving group 

Barrado y Navascues et al. 

1999 

Zuckerman et al. 200 

1), whose age has recently been de- 


bou ndary, isochrones with magnetic fields, and kinemat¬ 

ics (Binks fc Jeffries||2014 Malo et ak|2014 Mamajek &; 


Bell 


2014)7 

TT~ r Y7s) r 


IT 


F 27321 has a measured rota t ional velocity of ~ 120- 


130 kms 1 (Royer et al. 


2007 


Schroder et al. 


2009) 


and, while measurements of its metallicity are unreli¬ 
able, we can adopt the chemical composition of lower- 
temperature members of the same moving group. HIP 

10679 has a spectroscopic [Fe/H] = 0.07 ± 0.03, HIP 

10680 has [Fe/H] = 0. 09 ± 0.03, and HIP 2548 6 has 
[Fe/H] = 0.29 ± 0.03 (Valenti & Fischer| |2005[). As 
members of the same moving group, these stars should 
have nearly identical compositions. The inconsistency 
between their spectroscopic metallicities in the same sur¬ 
vey could indicate large systematic errors. The measure¬ 
ments do, however, hint at a slightly super-solar [Fe/H] 
for /3 Pic itself. 

The top panel of Figure [3] shows the age posterior prob¬ 
ability distribution for HIP 27321 under three metallicity 
priors. The red histogram, [Fe/H] = 0.1 ± 0.05, comes 
closest to the metallicities determined for later-type j3 
Pic members, and correctly indicates a young age for 
HIP 27321 itself. The lower limit on the star’s age is an 
artifact of our masking of the pre-main sequence. 

HIP 115738: HIP 115738 is an a 2 CVn variable star 
and a high- probability member of th e AB Doradus mov- 


Gagne et al. 

20141), with an age of ~100-150 Myr (Lull- 

man et al. 20 

05; Ortega et al. 2007 Barenfelcl et al. 2(T13). 

The metallicity of AB Dor is somewhat uncertain, with 


meas urements of [Fc/H] = 0.02 ± 0.02 ( Barcnfcld et al. 
M3| and [Fe/H] = 0.10±0.03 ([Biazzo et al.|2012| ), from 


spectroscopy of later-type members. 

Variable stars of a 2 CVn type are chemically peculiar 
with strong metal lines. This could imply a difference 
between the metallicity of the star and its atmosphere, 
an effect which we ignore. While HIP 115738’s variabil¬ 
ity is a generic problem for any isochrone analyses, the 
Tycho photometry is from a stack of 47 measurements 
(with a root-mean-square scatter in Vt of 0.029 mag). 
Stellar variability would be a much larger concern if we 
were to include non-simultaneous photometry at other 
wavelengths. 

HIP 115738’s variability does, in principle, permit a 
measurement of the its rotational period. Unfortunately, 
there seems to be little agreement on the star’s photo¬ 
metric penod_amongvariability surveys, ranging from 2 
hours ( Rimoldini et al.|2012 | to 1.4 days (Wraight et al. 


2012|. A rotation period of 2 hours is physically impossi¬ 


ble tor an A star of ~2 Rq. An undisputed measurement 
of the stellar rotation period would enable us to directly 
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constrain fl, the surface angular velocity. 

The second panel of Figure [3] shows the results of our 
analysis for HIP 115738. At a slightly super-Solar metal- 
licity, our posterior probability distribution is in excellent 
agreement with the age of the AB Dor moving group. 
HIP 115738’s age distribution would also skew somewhat 
younger if the star is a rapid rotator seen pole-on. 

HIP 13209: The late B star HIP 13209 (41 Ari) is 


a consensus member of AB Dor (Zuckerman et al.||2011 

Malo et al. 2013 

Gagne et al. 

2014). The star is a known 

spectroscopic bmarv and has 

been resolved at an angular 

separation of ~0'/l ( 

McAlister et al. 

1987 

). The binary 


5.2. The Pleiades Open Cluster 

The Pleiades open cluster is a large and extensively 
studied nearby stellar association. The Pleiades’ age is 
130 ± 20 Myr, as measured fr om the li thiu m depletion 
boundary (Barrado y Navascues et al.||2004|). It has also 
been estimate d at 100 Myr using st ellar isochrones with¬ 
out rotation (Meynet et al. 1993). Reddening for the 
bulk of the cluster is tl(B —V) = 0.04 mags, with a slight 
var iation from star to star (Vandenberg & Poll 1989f Tay¬ 
lor 20081. The distance to the Pleiades has been more 


controversial, but multiple lines of evidence, including 


and was not seen in near-infrared adaptive o ptics imaging 
(Roberts et al.l 2007 Janson et al. 2011). HIP 13209 


to ~135 pc (Pinsonneault et al. 

1998 Pan et al. 

2004 

Soderblom et al. 

2005 

Melis et al. 

2014|). As a 

large, 


et al. 2002). The metallicity of AB Dor is somewhat 
uncertain, with measu rements of [Fe/H] = 0.02 ± 0.02 
dBare nfeld et al.||2013 ) and [Fe/H] = 0.10 ±0.03 (Biazzo 
et al.||2012D, from spectroscopy for later-type members. 


has a measured rotational velocity of 175 kms (Abt test 0 f our statistical approach using the IGeorgv et al. 


HlP 13209 i s the earliest-type member of AB Dor listed 


Malo et al. (2013) and we include it here in spite of its 


companion. We make no attempt to remove the compan¬ 
ion’s flux, but note that any later-type companion would 
tend to bias our probability distributions to slightly older 
ages, i.e., make the star younger than our analysis im¬ 
plies. A G-type or later binary would contribute too little 
flux to make much difference. 

The third panel of Figure [3] shows the results for 
our three metallicity priors. With a Solar of slightly 
super-Solar metallicity, and neglecting the companion, 
our isochrone-based analysis agrees well with AB Dor’s 
known age. 

HIP 116805: The late-B star n And has been observed 
to host a substellar companion (Carson et al.|2013). The 
companion’s mass depends strongly on the host star’s 
age, which has been estimated from ~30 to ~300 Myr 
par son et al.||2013 Hinkley et al.|[2013 Bonnefoy et al. 


UI4p. The older age estimates "are - based on fitting 

nonrotating isochrones to the star’s position in color- 
luminosity space while the young age estimate relies on 
the sta r’s proposed membership in the Columba moving 
group (M alo et ah] [2013). We revisit the stellar parame¬ 
ters using the formalism we develop in this paper. 

We note that k And does have a spectroscopic metallic¬ 
ity m easurement giving [Fe/H] = —0.32±0.15 (Wu et al. 
2011). However, such a young, metal-poor star would 
be exceptional in the Solar neighborhood. Analyses of 
young clusters and moving groups, and of young stars 
in the local field, consistently find metallicities centered 
on or near the Solar value, with a dispersion ~0.1 dex 


(e.g. Gaidos & Gonzalez||2002 Biazzo et al.||2011 2012). 
We therefore use the same three metallicity priors as for 
the other stars discussed above: [Fe/H] = 0.1 ± 0.05, 
[Fe/H] = 0.0 ± 0.05, and [Fe/H] = -0.1 ± 0.05. The 
star has a rotational v sin i ps 160 kms -1 (Zorec & Royer 


2012 ). 

The bottom panel of Figure [3] shows our results for 
HIP 116805. A young age is excluded only if k And’s 
metallicity is strongly sub-Solar. If its metallicity is ap¬ 
proximately Solar or slightly super-Solar, HIP 116805’s 
position in color-magnitude space appears to be consis¬ 
tent with membership in Columba and a ^30 Myr age. 


(2013) isochrones. 

We restrict our analysis to late B Pleiads (~2.5 Mq), 
HIP 17527, 17588, 17664, 17776, 17862, and 17900. 
These are similar in mass to the stars that we will later 
use in analyses of the Hyades open cluster and the Ursa 
Majoris moving group. They have evolved sufficiently in 
^130 Myr to make an isochrone-based analysis possible, 
but are relatively far from the main-sequence turnoff. We 
adopt a parallax of 7.4 ± 0.2 mas (135 ± 4 pc) and a red¬ 
dening of E(B — V) = 0.04 mag, with Ry = 3.1, for 
all stars. We use a metallicity of [Fe/H] = 0.03 ± 0.05, 
as measur ed spectroscopically fo r 20 roughly Solar-type 
members (Soderblom et al. |2009| . 

Figure H shows our results) With a cluster [Fe/H] = 
0.03±0.057 the marginalized probability distribution is ~ 
95±35 Myr (ler), consistent with age determinat io ns us¬ 
ing th e lithi um depletion boundary (Barrado y Navascues 
et al. 2004) or higher-mass turnoff stars (Meynet et al. 
r993|). Fixing the metallicity to [Fe/H] = 0.03 has little 
effect on the results, changing the ler uncertainty from 
~35 to ^30 Myr. 

The most massive Pleiads, HIP 17499 (Electra) and 
HIP 17702 (Alcyone), both strongly prefer isochrone- 
based ages under 100 Myr (assuming the same distance 
and extinction), with Alcyone fitting ^6.5 Mq models at 
70 Myr. Electra’s best-fit age of ~90 Myr is con sistent 


with the other Pleiads. The Ekstrom et al. (20121 stellar 
models at D 0 /fl cr it = 0.58 predictTno stars A5 Mq near 
the main sequence: either Alcyone is a blue straggler, or 
the rotating stellar models fail to match massive B stars. 
Despite this caveat, the good agreement of the fit to six 
late B Pleiads with the cluster’s known age lends credi¬ 
bility to our application of the same ~2-2.5 Mq models 
to the Ursa Majoris and Hyades clusters. 

5.3. The Ursa Majoris Moving Group 

The Ursa Majoris moving group (UMa) is a coeval stel¬ 
lar association with an age o f 500 ± 100 Myr derived from 
isochrone fits (King et al.|2003). As a coeval association, 
it is also expected to be chemically homogeneous. We 
can therefore apply the fitting techniques presented in 
this paper to early-type UMa stars, jointly constraining 
the association’s age and metallicity. For this exercise, 
we choose five early-type sta rs that are listed as UMa 
members in King et al. (20031: HIP 53910, 59774, 6 2956, 
65477 , and 76267 We exclude HIP 6061, which |King| 
et al. (2003) list as a member but which lies well outside 





















































































































Fig. 4.— Posterior probability distribution in age and metallicity 
for the Pleiades, measured by fitting six late B stars, HIP 17527, 
17588, 17664 . 17776 . 17862, and 17900, to the stellar models of 
|Georgy et al.j (|2013|) with (red curves) and without (blue curves) 
allowing tor a stellar rotation. The cluster’s age has been mea- 
sured to be 130 it 20 Myr using the lithium depletion boundary 
(Barrado y Navascues e t ah 2004); this interval is shaded light 
blue. rnejGcorgy ct ai. (2UHs]l give a slightly younger age, albeit 
still consistent at Ict. There is little difference between the age and 


metallicity derived with and without considering rotation. As Fig¬ 
ure [2] shows, the rotating and nonrotating tracks for ~2-2.5 Mg 
stars only separate appreciably at older ages. 


the cluster core and is inconsistent with the other stars’ 
ages. In the case of HIP 76267, we subtract the flux of its 


eclipsing binary companion (a G5V star, Tomkin fe Pop- 
per|l986|. HIP 65477 has an M-type companion (|Mama- 


jek et ai. 


2010) which is too faint to appreciably affect 
its optical flu x. We exclude the other early-type UMa 
members from King et al. (2003), all of which appear to 
be spectroscopic binaries with companions of unknown 
spectral type. Four of our five stars, HIP 53910, 59774, 
62956, and 65477, are members of the UMa nucleus. 

When fitting the ensemble of stars, we assume a Gaus¬ 
sian metallicity prior for the cluster centered on [Fe/H] = 
0 with a dispersion of 0.1 dex. We assume the individ¬ 
ual stars to have the (unknown) cluster value. Each in¬ 
dividual star may have its own parallax, rotation and 
orientation; we marginalize over all of these parameters 
star-by-star. We then multiply the joint posterior prob¬ 
ability distributions in age and metallicity for all stars, 
obtaining the results shown in Figure [5j The upper-left 
density plot assume all stars to be nonrotating while the 
lower-right assumes a Maxwellian prior on f2o/n cr i t . The 
inner and outer contours enclose 68% and 95% of the 
probability, respectively. The dot-dashed curves show 
the results if the metallicity prior is a delta function at 
[Fe/H] = 0.03. 

When using the rotating models, all of the stars are 
consistent with one another in their age-metallicity con¬ 
straints; their la contours overlap one another. The 
derived age is in good agreement with the 500 ± 100 


Fig. 5.— Joint posterior probability distribution for the age and 
metallicity of five early-type Ursa Majoris moving group members, 
with a metallicity prior of [Fe/H] = 0 it 0.1. Our derived age wi th 
rotation agrees well with ~500 estimated by |King et al. 2003b. 
The red dot-dashed curves show the constraints from assuming a 
delta-function metallicity prior [Fe/H] = 0.03, the central value 
meas ured from FGK candidate UMa members (pTabernero et al. 
2014). When neglecting rotation, the age/metallicity constraints 
tor HIP 62956 do not agree with the other stars; multiplying the 
probability densities yields the implausibly young age and high 
metallicity shown. 


Myr reported by King et al. 


(2003). Our isochrone anal¬ 


ysis also favors a very slightly super-Solar metallicity, 
in agreement with rec ent spectroscopic measur ements of 


later-type members. Tabernero et al. (2014), using a 


sample of 44 candidate FGK members, find that 29 of 
the 44 stars share a similar chemical composition, with 
[Fe/H] = 0.03 ± 0.07 dex. Using a delta function metal¬ 
licity prior at [Fe/H] = 0.03 (dot-dashed red lines) gives 
an age of 530 ± 40 Myr (2er formal errors). 

The agreement is much worse when neglecting rota¬ 
tion. Much of this is due to HIP 62956: its ler con¬ 
tour does not overlap those of the other stars in age- 
metallicity space. Excluding this star yields good agree¬ 
ment in the nonrotating case at an age of just over 400 
Myr, which is still ~20% younger than in the rotating 
case. Including or excluding HIP 62956 has a negligible 
effect on the constraints from rotating stellar models. 

5.4. The Hyades Open Cluster 

We now revisit the age of th e Hyades open cluster cal¬ 
culated to be 625 ±50 Myr by [Perryman et al.|(1998| us¬ 
ing nonrotating isochrones with convective overshooting. 
The|Perryman et al. (19981 constraint is based mostly on 
five stars near the mam-sequence turnoff without indica¬ 
tions of multiplicity: HIP 20542, 20635, 21029, 21683, 
and 23497. Reddening be tween the Sun and the Hyades 
is negligible (Taylor 2006). We use the same five stars in 
our present framework to place new constraints on the 
Hyades age, adopt ing a metallicity of [Fe/H] = 0.10 (Tay- 
lor fe Joner||2005 | with a conservative Gaussian error of 
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0.05 dex. We also add the individual components of the 
binary HIP 20894 as two additional stars, with a magni¬ 
tude difference AVr = 1.10 and A (Bt — Vt) = —0.006 
(Peterson et al. 1993). As for Ursa Majoris, we assume 
the stars to share a composition and marginalize over 
mass, rotation, and orientation separately for each star. 

The lengthening of the main-sequence lifetime with ro¬ 
tation maps onto an older inferred age of ~750 Myr for 
the Hyades, shown in Figure [6] As Figure [2] shows, a ro¬ 
tating isochrone can also overlap a slightly more massive 
nonrotating isochrone near the main sequence turnoff. 
Indeed, the difference between the maximum likelihood 
masses with and without including rotation is typically 
~0.05-0.1 Mq, which corresponds to a ~10% difference 
in main-sequence lifespan. These effects combine to in¬ 
crease the Hyades age from just under 600 Myr to about 
750 Myr. We note that the nonrotating models used by 


Perryman et al. (19981 inclu ded slightly s t ronge r core 
convective overshoot than the Georgy et al. (20131 mod¬ 
els (a = 0.2 compared to a = 0.1 j. This mimics some 
of the effects of rotation, and accounts for the modest 
difference between the nonrotating models here and in 


Perryman et al. (1998). 

All seven turnoff stars, including the two components 
of HIP 20894, are consistent with an age of ~750-800 
Myr. The 800 Myr, Qo/f2 cr jt = 0.58 Zq isochrone in 
Ekstrom et al. (2012) has its red clump at B — V ss 1, 
and V ss 0.4, completely consistent with the location 
of the four r ed giants in Fi g ure 2 1 of | Perry man et ah] 
(|1998|). The |Georgy et al.| (|2013D nonrotating models 
strongly prefer a high metailicity for the cluster, signifi¬ 
cantly higher than the spectroscopic value. The rotating 
models are fully consistent with the spectroscopic metai¬ 
licity; fixing Fe/H] =0.1 gives an age of 800 ± 50 Myr. 

Open clusters like the Hyades are often used as bench¬ 
marks to calibrate age dating techniques for lower mass 
stars. If the cluster’s age really is closer to 800 Myr than 
to 600 Myr, it may require a re-calibration of seco ndary 


age indicators like st ellar rotation and activity (Mamajek 
fc Hillenbrand|[2008 ). At the very least, our results sug¬ 


gest that the uncertainty in the Hyades’ age is underes¬ 
timated. Resolving the discrepancy might require other 
age indicators for these stars or additional tests of the 
rotating stellar models. Asteroseismology of lower-mass 
stars provides a relatively direct measurement of the cen¬ 
tral sound speed, and by extension, the mean molecular 
weight and helium fraction. The T ransiting Exoplanet 
Survey Satellite (TESS, |Ricker et al.||2014|) could soon 
provide crucial data for F and G-type Hyades members. 

Two of our seven turnoff stars, HIP 20894 A and B, 
form a binary with dynamically measured masses. The 
most recent and most precise masses are 2.86 ± 0.06 Mq 
and 2.16 ± 0.02 Mq for the p rimary and secondary, re¬ 
spectively (Torres et al.pOll). If these masses are cor¬ 
rect, they raise more questions than they answer: none 
of our 2.86 Mq stellar models come within 1 mag of the 
observed lum inos ity of HIP 20894 A at any age. Fig¬ 
ure 10 of Torres et al. ([2011) confirms this, and further 
shows that 2.86 Mq stars with HIP 20894 A’s colors 
are rapidly traversing the subgiant branch. We note 
that our analysis, with a ~2.4 or 2.5 Mq HIP 20894 
A, places the star on the main-sequence turnoff before 
the subgiant branch, where it is evolving ten times more 
slowly through color-magnitude space. Previous orbital 



Age (Myr) 


Fig. 6.— Age-metallicity constraints on the Hyades open cluster 
using seven main-sequence turnoff members (HIP 20542, 20635, 
21029, 21683, 23497 and the two components of HIP 20894). The 
inner and outer contours enclose 68% and 95% of the probability, 
respectively; the left contours and blue curves neglect rotation. 
The best-fit age is ~750 Myr when including an uncertainty of 
0.05 in [Fe/H] (solid red line), or ~800 Myr with at [Fe/H] = 0.1 
(dot-dashed red curve). These ages are significantly older than the 
consensus age of 625 Myr , which was calculated largely based on 
five of these seven stars ([Perryman et al.||1998). We shade the 
previous age constraint in light blue tor comparison. 


solutions of the same system have found very d iffe rent 


masses, albeit with m uch less data (Torres et al. 


Armstrong et al. 12006 1. A confirmation of HIP 20894 A’s 


1997 


high mass would present major challenges to models of 
stellar evolution. 


6. CONCLUSIONS 

In this paper, we have used recently computed stellar 
m odels with an imp roved treatment of stellar rotation 
( Georgy et al.|2013 ) to derive stellar parameters, partic¬ 
ularly ages, of nearby earl y-type stars. We combin e AT- 
LAS9 model atmospheres ( iCastelli fe Kurucz|2004| with 
the two-dimensional model of l' e s and log g m LR11 to 
compute synthetic magnitudes as a function of orienta- 
tion. We then use fine grids of nonrotating isochrones 
(Girardi et al. 2002) to interpolate between the rotat¬ 
ing models in time, mass, and metailicity. In the case of 
metailicity, we also use the nonrotating models to extrap¬ 
olate the rotating models to super-Solar metallicities. 

We have applied a Bayesian age analysis to obtain 
posterior probability distributions of the stellar param¬ 
eters assuming random orientations and a standard ini¬ 
tial mass function. While we have computed synthetic 
photometry through a wide range of filters, the analy¬ 
sis we present relies only on the Tyc ho BtVt system 
and Hipparcos parallaxes (ESA 1997). The phot ome¬ 
try and astrometry, which have s ince been updated (Hpg 
et al. 2000 van Leeuwen 2007), provide uniform and 
high-quality data tor ail bright stars. 

We have applied our analysis first to stars with known 
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ages, including the Pleiades open cluster. Our use of the 
new rotating isochrones recovers, within its errors, the 
known ages of /3 Pictoris and (with caveats) two mem¬ 
bers of the AB Doradus moving group. Fitting of the late 
B Pleiads recovers an age of ~100 Myr, in good agree¬ 
ment with previous results using both i soch r one fitting 
arid tire lithium d epleti o n bo undary (Meynet et al.|1993 


Barrado y Navascues et al. 2004). Fitting the massive 
Pleiad Alcyone favors a younger age. This could indi¬ 
cate that Alcyone is a blue straggler, or that the rotating 
stellar models fail at higher masses. 

Fitting ~2-2.5 M 0 models to the Ursa Majoris moving 
group agrees well wit h previously published ages for that 
cluster)King et ah [2003). Fitting the turnoff stars of the 
Hyades, however, implies a cluster age of ~750-800 Myr, 
significantly older than the previous estimate of 6 25 ± 50 
Myr using the same stars (Perryman et al. 1998). This 
difference results from an increase in both mam-sequence 
lifetime and luminosity due to rotational mixing. 

The older age we derive for the Hyades turnoff stars 
raises questions about either the cluster’s age, stellar 
modeling, or both. Including stellar modeling uncertain¬ 
ties could imply a much larger error of at least 100 Myr 
for the cluster age, which would propagate into secondary 
age indicators calibrated on clusters like the Hyades. Re¬ 
solving these modeling uncertainties will likely require 
independent age indicators, perhaps asteroseismology of 
F and G type Hyades stars. 


Finally, we provide a web interface at 
www.bayesianstellarparaLmeters.info where a 
user may obtain the posterior probability distributions 
of age, mass, and orientation for any bright ( H P < 9), 
blue (B — V< 0.25) Hipparcos star with a good distance 
( voIo^j > 15). The web server currently assumes negli¬ 
gible extinction and uses only Tycho BtVt photometry 
and Hipparcos parallaxes. A careful analysis would 
be needed to place other photometry, like 2MASS, on 
the same absolute scale, and ensure that the resulting 
ages remain consistent with one another and with 
independently determined ages. 


The authors would like to thank the computing staff at 
the Institute for Advanced Study for help in setting up 
a web interface, Waqas Blratti for the server script tem¬ 
plate, and Josh Schlieder, Lynne Hillenbrand, and Geza 
Kovacs for helpful comments on the manuscript. The 
authors thank the referee, Juan Zorec, for helpful com¬ 
ments and for pointing out an error in the age scale. This 
work was performed in part under contract with the Jet 
Propulsion Laboratory (JPL) funded by NASA through 
the Sagan Fellowship Program executed by the NASA 
Exoplanet Science Institute. This research has made use 
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REFERENCES 


Abt, H. A., Levato, H., & Grosso, M. 2002, ApJ, 573, 359 
Anders, E., Sz Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 
197 

Armstrong, J. T., Mozurkewich, D., Hajian, A. R., et al. 2006, 
AJ, 131, 2643 

Asplund, M., Grevesse, N., Sauval, A. J., Sz Scott, P. 2009, 
ARA&A, 47, 481 

Aufdenberg, J. P., Merand, A., Coude du Foresto, V., et al. 2006, 
ApJ, 645, 664 

Barenfeld, S. A., Bubar, E. J., Mamajek, E. E., Sz Young, P. A. 

2013, ApJ, 766, 6 
Barnes, S. A. 2003, ApJ, 586, 464 
Barnes, S. A. 2007, ApJ, 669, 1167 

Barrado y Navascues, D., Stauffer, J. R., Sz Jayawardhana, R. 
2004, ApJ, 614, 386 

Barrado y Navascues, D., Stauffer, J. R., Song, I., Sz Caillault, 
J.-P. 1999, ApJ, 520, L123 
Bessell, M., & Murphy, S. 2012, PASP, 124, 140 
Biazzo, K., D’Orazi, V., Desidera, S., et al. 2012, MNRAS, 427, 
2905 

Biazzo, K., Randich, S., Palla, F., Sz Briceno, C. 2011, A&A, 530, 
A19 

Binks, A. S., &; Jeffries, R. D. 2014, MNRAS, 438, Lll 
Bonnefoy, M., Currie, T., Marleau, G.-D., et al. 2014, A&A, 562, 
Alll 

Carson, J., Thalmann, C., Janson, M., et al. 2013, ApJ, 763, L32 
Casagrande, L., Schonrich, R., Asplund, M., et al. 2011, A&A, 
530, A138 

Castelli, F., Sz Kurucz, R. L. 2004, ArXiv Astrophysics e-prints, 
astro-ph /0405087 
Chabrier, G. 2003, PASP, 115, 763 

Chiappini, C., Romano, D., & Matteucci, F. 2003, MNRAS, 339, 
63 

Claret, A. 1998, A&AS, 131, 395 

Connelly, J. N., Bizzarro, M., Krot, A. N., et al. 2012, Science, 
338, 651 

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, 2MASS 
All Sky Catalog of point sources. 

Dotter, A., Chaboyer, B., Jevremovic, D., et al. 2008, ApJS, 178, 
89 

Ekstrom, S., Georgy, C., Eggenberger, P., et al. 2012, ASzA, 537, 
A146 

ESA, ed. 1997, ESA Special Publication, Vol. 1200, The 
HIPPARCOS and TYCHO catalogues. Astrometric and 
photometric star catalogues derived from the ESA 
HIPPARCOS Space Astrometry Mission 


Espinosa Lara, F., Sz Rieutord, M. 2011, A&A, 533, A43 
Gagne, J., Lafreniere, D., Doyon, R., Malo, L., Sz Artigau, E. 
2014, ApJ, 783, 121 

Gaidos, E. J., Sz Gonzalez, G. 2002, New Astron., 7, 211 
Georgy, C., Ekstrom, S., Granada, A., et al. 2013, A&A, 553, A24 
Georgy, C., Granada, A., Ekstrom, S., et al. 2014, ASzA , 566, A21 
Girardi, L., Bertelli, G., Bressan, A., et al. 2002, A&A, 391, 195 
Glatzmaier, G. A. 1985, ApJ, 291, 300 

Hempelmann, A., Schmitt, J. H. M. M., Schultz, M., Ruediger, 

G., Sz Stepien, K. 1995, A&A, 294, 515 
Hinkley, S., Pueyo, L., Faherty, J. K., et al. 2013, ApJ, 779, 153 
Hpg, E., Fabricius, C., Makarov, V. V., et al. 2000, ASzA, 355, L27 
Janson, M., Bonavita, M., Klahr, H., et al. 2011, ApJ, 736, 89 
King, J. R., Villarreal, A. R., Soderblom, D. R., Gulliver, A. F., 

& Adelman, S. J. 2003, AJ, 125, 1980 
Kroupa, P. 2001, MNRAS, 322, 231 

Lagrange, A.-M., Gratadour, D., Chauvin, G., et al. 2009, A&A, 
493, L21 

Lebreton, Y., Goupil, M.-J., Sz Montalban, J. 2014, ArXiv 
e-prints, 1410.5336 

Li, C., de Grijs, R., Sz Deng, L. 2014, Nature, 516, 367 
Luhman, K. L., Stauffer, J. R., Sz Mamajek, E. E. 2005, ApJ, 

628, L69 

Malo, L., Doyon, R., Feiden, G. A., et al. 2014, ApJ, 792, 37 
Malo, L., Doyon, R., Lafreniere, D., et al. 2013, ApJ, 762, 88 
Mamajek, E. E., Sz Bell, C. P. M. 2014, ArXiv e-prints, 1409.2737 
Mamajek, E. E., Sz Hillenbrand, L. A. 2008, ApJ, 687, 1264 
Mamajek, E. E., Kenworthy, M. A., Hinz, P. M., & Meyer, M. R. 
2010, AJ, 139, 919 

McAlister, H. A., Hartkopf, W. I., Hutter, D. J., Sz Franz, O. G. 
1987, AJ, 93, 688 

Melis, C., Reid, M. J., Mioduszewski, A. J., Stauffer, J. R., Sz 
Bower, G. C. 2014, Science, 345, 1029 
Meynet, G., Sz Maeder, A. 2000, A&A, 361, 101 
Meynet, G., Mermilliod, J.-C., Sz Maeder, A. 1993, A&AS, 98, 477 
Nielsen, E. L., Liu, M. C., Wahhaj, Z., et al. 2013, ApJ, 776, 4 
Nieva, M.-F., Sz Przybilla, N. 2012, A&A, 539, A143 
Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., 
&; Vaughan, A. H. 1984, ApJ, 279, 763 
Ortega, V. G., Jilinski, E., de La Reza, R., Sz Bazzanella, B. 2007, 
MNRAS, 377, 441 

Pan, X., Shao, M., Sz Kulkarni, S. R. 2004, Nature, 427, 326 
Parker, E. N. 1955, ApJ, 122, 293 

Perryman, M. A. C., Brown, A. G. A., Lebreton, Y., et al. 1998, 
A&A, 331, 81 










11 


Peterson, D. M., Stefanik, R. P., & Latham, D. W. 1993, AJ, 105, 
2260 

Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 
612, 168 

Pinsonneault, M. H., Stauffer, J., Soderblom, D. R., King, J. R., 
& Hanson, R. B. 1998, ApJ, 504, 170 
Preston, G. W. 1974, ARA&A, 12, 257 

Przybilla, N., Nieva, M.-F., & Butler, K. 2008, ApJ, 688, L103 
Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, in Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 9143, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, 20 
Rimoldini, L., Dubath, P., Siiveges, M., et al. 2012, MNRAS, 427, 
2917 

Roberts, L. C., Jr., Turner, N. H., & ten Brummelaar, T. A. 

2007, AJ, 133, 545 

Roy, J.-R., & Kunth, D. 1995, A&A, 294, 432 

Royer, F., Zorec, J., & Gomez, A. E. 2007, A&A, 463, 671 

Salpeter, E. E. 1955, ApJ, 121, 161 

Schroder. C., Reiners, A., & Schmitt, J. H. M. M. 2009, A&A, 

493 1099 

Sestito, P., & Randich, S. 2005, A&A, 442, 615 
Smith, K. C. 1996, Ap&SS, 237, 77 
Soderblom, D. R. 2010, ARA&A, 48, 581 

Soderblom, D. R., Laskar, T., Valenti, J. A., Stauffer, J. R., & 
Rebuff, L. M. 2009, AJ, 138, 1292 
Soderblom, D. R., Nelan, E., Benedict, G. F., et al. 2005, AJ, 

129, 1616 

Tabernero, H. M., Montes, D., Gonzalez Hernandez, J. I., & 
Ammler-von Eiff, M. 2014, ArXiv e-prints, 1409.2348 


Takeda, G., Ford, E. B., Sills, A., et al. 2007, ApJS, 168, 297 
Taylor, B. J. 2006, AJ. 132, 2453 
Taylor, B. J. 2008, AJ, 136, 1388 
Taylor, B. J., & Joner, M. D. 2005, ApJS, 159, 100 
Tomkin, J., & Popper, D. M. 1986, AJ, 91, 1428 
Torres, G., Stefanik, R. P., & Latham, D. W. 1997, ApJ, 485, 167 
Torres, K. B. V., Lampens, P., Fremat, Y., et al. 2011, A&A, 525, 
A50 

Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 
van Belle, G. T., Ciardi, D. R., ten Brummelaar, T., et al. 2006, 
ApJ, 637, 494 

van Leeuwen, F. 2007, A&A, 474, 653 
Vandenberg, D. A., & Poll, H. E. 1989, AJ, 98, 1451 
von Zeipel, H. 1924, MNRAS, 84, 684 
Wilson, O. C. 1963, ApJ, 138, 832 

Wraight, K. T., Fossati, L., Netopil, M., et al. 2012, MNRAS, 

420, 757 

Wu, Y., Singh, H. P., Prugniel, P., Gupta, R., & Koleva, M. 2011, 
A&A 525 A71 

Yi. S., Demarque, P., Kim, Y.-C., et al. 2001, ApJS, 136, 417 
Za'hn, J.-P. 1992, A&A, 265, 115 
Zorec, J., & Royer, F. 2012, A&A, 537, A120 
Zuckerman, B., Rhee, J. H., Song, I., & Bessell, M. S. 2011, ApJ, 
732, 61 

Zuckerman, B., Song, I., Bessell, M. S., & Webb, R. A. 2001, 

ApJ, 562, L87 



